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Abstract 

Isospin-breaking in B — > tttt caused by ir° — n — rf mixing is studied in a model-independent way 
using flavor SU(3). Measured branching ratios for B + — ► 7r + 7r°, B + — ► 7r + ?/') and B° — > TT°n^'^ 
imply an uncertainty in a smaller than 1.4°. We find a negligible effect of ir° — n — rf mixing 
on a in B — > p7r. Characterizing the effect of p° — a; mixing in B — > pp and in B — > p-zr by the 
two-pion invariant mass dependence, we point out a way of constraining this effect experimentally 
or eliminating it altogether. We show that a model- independent shift in a caused by electroweak 
penguin amplitudes in B — > 7T7r and B — ► pp, Aoewp = (1-5 ± 0.3)°, may be slightly different in 
B — > p7r. Other sources of isospin-breaking in these processes are briefly discussed. 
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I. INTRODUCTIION 



Isospin symmetry provides triangle relations for B — > tttt and B — > ttit, which are governed 
by / = and 1 = 2 amplitudes, 

+ V2A 00 - V2A +0 = , Aij = A(B ^ ttV) , (1) 
A + _ + V2A 00 - >/2A_„ = , A^ = A(B° -> ttV) . (2) 

These relations enable an extraction of the phase a = 2 ee Arg(— V t * b V t d/V* b V ud ) from time- 
dependent CP asymmetries, SVtt and C n7T , in B°(t) — > P. The asymmetries determine 
sin2a eff , 

sin(2a e ff) = — , . (3) 

The shift Aa ee a cff — a caused by the penguin amplitude is given by 

Aa ee ^Arg (e 2i ^A + _A* + _) . (4) 

We define a measurable phase A«o, given in terms of angles in the B and B triangles, 
<j) ee Arg (A + _A* +0 ), ee Arg {A + _A%): 

Aa ee 1(0 - 0) = i [Arg (e 2 ^ + _A;_) - Arg (e 2 ^A_ A; o )] . (5) 

Neglecting very small electroweak penguin amplitudes which will be discussed below, a phase 
relation holds between the two AI = 3/2 tree amplitudes, 

A +0 = e 2i ^A_ . (6) 

This implies Aa = Acto, fixing the relative orientation of the B triangle and the B triangle 
(rotated by 2 / y) such that these sides overlap. In this configuration Aa is half the angle 
between sides corresponding to B° — ► 7r + 7r~ and B° — > 7r + 7r~. This determines Aa up 
to a fourfold ambiguity (including a sign ambiguity) related to the four possible relative 
orientations of the two triangles. In case that \A(B° — > 7r°7r°)| and \A(B° — > 7r°7r°)| are 
not separately measured, while the charge-averaged neutral pion rate is measured, one may 
obtain upper bounds on |Aa| Q. 

The same method applies to polarization states in B°(t) — > p + p~ which are CP- 
eigenstates, in particular to an even-CP longitudinally polarized state which is found to 
dominate this process [3j. A variant of isospin symmetry can also be used to learn a in 

and 



dominate tms process [o|j. A variant ol isospin symmetry can also be used to learn . 
B°(t) — > pn Application of these methods to recent measurements by BaBar 0,0] 

Belle [7], where B — ^p^played a dominant role, provides the currently most accurate direct 

by resolving 



>P£Piayed a dominant role, provides tne currently most ac 
determination of a jaM, ol = (100l^ )°. This precision can be improved 1(J 
the sign ambiguity in Aa under reasonably mild and testable assumptions about strong 
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phase differences between tree and penguin amplitudes in these processes. The error in a is 
expected to be reduced further by improving the measurement of_the direct CP asymmetry 
in B° -> ttV 



11 1, or by using a prediction of a SCET analysis [3| that the phase difference 
between tree and color-suppressed amplitudes in B + — > 7r + 7r° is small 

At this level of precision one is required to consider small electroweak penguin (EWP) 
amplitudes and corrections from isospin breaking caused by the u and d charge and mass 
differences. These corrections modify the geometry of the B and B amplitude triangles. A 
model-independent study of electroweak penguin contributions in B — ► tttt was performed 
in fl4 1 . Instead of overlapping with each other, the two sides of the two triangles, A +Q and 
e 2i7 A_o, were shown to form a calculable relative angle. Neglecting EWP operators with 
small Wilson coefficients [c-j and eg), isospin symmetry relates dominant AI = 3/2 EWP 
operators to AI = 3/2 current-current operators in the effective Hamiltonian, implying 

(Aa - Aa„) EWP = \kr g {^ A^A* +0 ) = -§ f^±^) ^^sina 

A A \ C\ + C 2 / | V u b v ud\ 

3 f cg + ci \ sin(/3 + a) sin a _ ^ ^ 3 sin(/3 + a) sin a 



2 \ ci + C2 / sin /3 sin f3 

The measured values of j3 and a |8|, /3 = (23.3 ± 1.6)°, a = (100l^ )°, lead to a small 
calculable value with a negligible error, (Aa — A«o)ewp — (1-5 ± 0.3)°. This shift must be 
included in the determination of a by using a = a c s — Aa — (Aa — A«o)ewp- 

Isospin-breaking due to nonzero u and d quark mass and charge differences has several 
effects on the analysis of B — >• tttt, B —>■ pp and B — > pn. An important effect in B — ► Tin, 
caused by 7r° — 77 — 7/ mixing, was studied several years ago by Gardner [l^ using generalized 
factorization She concluded that the resulting error on the extracted value of a in the 
above range is about 5° including EWP contributions. The uncertainty may be even larger 
due to the approximation involved in this estimate. This would limit severely the future 
accuracy of determining a in B — > tttt. Gardner and Meissner 0| discussed briefly the 
appearance of a small AI = 5/2 amplitude in B — » tttt which violates the isospin triangle 
relation. They mentioned isospin violation in B — > pn caused by tt° — rj — rf mixing and 
by p — uj mixing, pointing out that the presence of an additional AI = 5/2 amplitude in 
B — > pir, involving the same weak phase as the tree amplitude, does not affect the isospin 
analysis. Refs. 18, 0, 0] studied direct CP violation in B + >° — ► (7r + 7r~) Pja; 7r + ' and in 



B + — > (tt + tt~) p ^p + caused by p — ui mixing. These CP asymmetries affect the analyses 
of isospin related processes. Ref. 0] studied numerically the uncertainty in determining a 
in B — > pp, assuming that isospin violating corrections in tree and penguin amplitudes in 
\/2A(B + — ► p + p°) are at a level of 4% relative to tree and penguin amplitudes in B° — > p + p~ ■ 
The purpose of this work is to analyze in a model-independent manner isospin-breaking 
effects, in particular the effects of tt° — rj — 77' mixing and p — uo mixing, on determining a 
in B — > tttt, B — > pp and — > p7r. In Section II we will apply flavor SU(3) to 5 decays 
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into two charmless pseudoscalars, relating isospin-breaking terms in B — > txtx to amplitudes 
of B — > 7T7] and B — > -nrf . Using measured rates, we will show that the effect of tt° — rj — rf 
mixing on determining a in B — > im is considerably smaller than estimated by Gardner. 
Turning in Section III to discuss the effects of p — to mixing on determining a in B — > pp, we 
will show how to include this effect experimentally, without having to rely on a calculation 
of p — u mixing parameters. Section IV studies 7r° — rj — rf mixing, p — uj mixing and the 
effect of electroweak penguin amplitudes in B —>■ p-n . In Section V we discuss briefly other 
sources for isospin-breaking, while Section VI concludes. Appendix A presents experimental 
constraints on parameters describing p — uj mixing. 



II. EFFECT OF vr° - rj - rf MIXING IN B -> vrvr 

The mixing of tt° with rj and rf introduces isospin-breaking in B — ► ttti through an 
additional 1 = 1 amplitude, while the isospin conserving terms obey the triangle relation 
(UJ). We will use flavor SU(3) symmetry to estimate the isospin-breaking terms. SU(3) 
breaking corrections and smaller annihilation-like amplitudes which we neglect in these 
terms are higher order, and are expected to introduce an uncertainty at a level of 30%. 
A convenient way of applying flavor SU(3) to charmless B decays into two pseudoscalars 
is in terms of graphical representations describing flavor flow topologies 2l|, |22|, |23|, \24\. 
SU(3) amplitudes for two octets in the final state consists of a "tree" amplitude (t) a "color- 
suppressed" amplitude (c) a "penguin" amplitude (v). Three annihilation-like amplitudes 
(a, e and pa) are expected to be much smaller 12 3L 1251 ] and will be neglected. The remaining 
three amplitudes contain EWP contributions |26j], the overall effect of which can be taken 
into account as summarized in Eq. ((Zj). This effect can be included as explained above, and 
will therefore be disregarded in this section. For a singlet and an octet in the final state 
one has three SU(3) amplitudes j2l|, of which a "singlet penguin" amplitude (s) dominates, 
while two annihilation-like amplitudes will be neglected [24| . 



We use quark content for mesons as in [2J,|2J|, but a somewhat different phase convention: 

B° = —db, B + = —ub, 7r 3 = —=(uu — dd), tt + = —ud, it~ = du, 

V2 

1 _ 1 

rj = —=(uu + dd — ss), rj' = —=(uu + dd + 2ss) . (8) 
V3 v6 

The rj and rf correspond to octet-singlet mixtures, 

rj = rjs cos 9 — rji sin 9 , rf = r/ 8 sin 9 + r\\ cos 9 , 

r] 8 = —={uu + dd — 2ss) , r)i = —=(uu + dd + ss) , (9) 
V6 V3 

with an "ideal" mixing angl e 9 = 9q = sin _1 (— 1/3) = —19.5°. A slightly larger magnitude, 
9 = —22°, was used in [l5(, while a slightly smaller magnitude, 9 = (—15.7 ± 1.7)°, was 
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obtained in a very recent phenomenological fit [2^. While in the most part we use the value 
6*0, at the end of this section we discuss briefly the effect of a variation in 9. 

Expressions for decay amplitudes in terms of graphical SU(3) contributions, for final 



states involving pairs of isotriplet pions, and 
rj', are obtained in a straightforward manner 



or pairs involving an isotriplet pion and rj or 



23, 



24, 



A 3v = 



-t+p 
1 

1 

7! 



^33 

(2p + s) , 
(t + c + 2p + s 



1 



T2 [t + C) 



A 



3r?' 



v/3 



(p + 2s) , 



A, 



(t + c + 2p + 4s) . 



(10) 



The first three amplitudes for pure isotriplet pions obey clearly the isospin triangle relation 
(JTJ). The purely AJ = 3/2 amplitude A +3 = (t + c)/y/2 has a weak phase 7. 

The mixing of 7r°, r] and 7/ introduces a small isospin singlet component into the domi- 
nantly isotriplet neutral pion state, 

|7r ) = |7r 3 )+ e |r / ) + e '|V) . (11) 

Values e = 0.014, e' = 0.0077 were used by Gardner [l^, based on a calculation apply- 



ing chiral perturbation theory |29j| . We will take ranges of values as obtained in a recent 
update 0, e = 0.017 ± 0.003, e' = 0.004 ± 0.001. 

Neglecting terms quadratic in e and e', we find decay amplitudes for the neutral pion 
state given by (JTTJ): 



.4 



+0 



A 



00 



A +3 + eA +v + e'A +rf 

1 1 \p2* 

— (t + c)(l + e ) + -^=£(2^ + s) + —e\p + 2s) 

A 33 + v^eAg,, + v^e'Aav 



! 1 \j2i 

■y=(c - p) + -^e(2p + s) + -j=e'{p + 2s) 



where 



e 



e + 



1 £ / 



0.016 ±0.003 . 



(12) 



(13) 



(14) 



We note the factors a/2 in the first line of Eq. (|13|). This takes into account final states of 
identical particles in A 00 and A 33 [31], compared to states which must be symmetrized in 
A 3r] and A 3r) >. Squares of amplitudes give decay rates when common phase space factors are 
implied. 

Using these expressions, we find two important consequences of this amplitude decom- 
position which includes tt° — rj — rf mixing: 
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1. The triangle relation (0) is modified only slightly: 

+ V2A 00 - V2A +Q {1 - e ) = 



(15) 



2. The amplitude A + q can be written in terms of the pure AI = 3/2 amplitude, A + %, 
carrying a weak phase 7, corrected by isospin-breaking terms involving A 0r) and A 0r) i, 



A +0 = A +3 (l + eo) + v^eA^ + v^e'A 



(16) 



Our first conclusion is therefore that the physical B — > ttit and B — > 7T7t decay amplitudes 
still obey triangle relations. The isospin-breaking factor, 1 — eo, multiplying the amplitude 
A +0 in (fT5]). can be absorbed in this measurement. Since e is calculated to be between 
one and two percent, while the current error in \A + q\ is about 5% (see Eq. (|2())l below), the 
factor 1 — eo starts to play a non-negligible role and must be included in the construction of 
the isospin triangles |22|. The remaining error from the theoretical uncertainty in eo given 
in (JT3j) is only a fraction of a percent, causing a negligible error in determining Acto from 
the angles in the two isospin triangles. 

The second result, Eq. (fTfij). implies that A +0 and its charge-conjugate no longer obey 
the exact phase relation (JHJ). That is, since the weak phases of the small isospin-breaking 
terms in ()16j) differ from 7, the triangle (fTSJl and the corresponding triangle for B amplitudes 
rotated by an angle 27 do not share exactly a common base, A + o 7^ e 2ry A_o- Denoting 



Arg 



A, 



A, 



Are 



A 



A- 



(17) 



this introduces a change, Act — Aao, given to first order in e and e' by 
(Aq - Aao).-,- V = ^Arg(e 2l ^_o^; o ) 



(18) 



1 



V2\A, 



el \A 0r] \ smifj v — \A 0v \ smi/) v ) + e'f \A 0ri i\ sin^/ — {A^l smip n 



Given that the phases if) o and ^ o are unknown, an immediate upper bound on \Aa 
Aa \ may be obtained by taking if) o = —ip„c) = 7r/2: 



|(Aa - Aa ) n . 



< e 



< 



\Aq v 


+ \A r,\ 


V2\A +0 \ 




\ /SO, 



+ e' 



\Ao v ' 


+ 


A 


Or?' 


V2\A. 


fO 





T 



B 



+ e 



+0 



(19) 



Here Bij = (|^4ij| 2 + \Aij\ 2 )T B /2 denote charge-averaged branching ratios for corresponding 
decays, and t + /tq is the lifetime ratio of B + and B°. We neglect tiny corrections (at a level 
of a percent) in phase space factors. 
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Using world averaged values |33l |. 
7b 



1.081 ± 0.015 , B +0 = (5.5 ± 0.6) x 10" b [3J, |3J , 

B 0r , < 2.5 x 10" 6 (90% CL) 0,0 , Bty < 3-7 x 10^ 6 (90% CL) , (20) 



we find at 90% CL 



| (Ac* - Aao)^-^-^' I < 105e + 1.28e' = 1.6° . 



(21) 



The phases ip v and if) v i may actually be measured within discrete ambiguities through 
two triangle relations implied by Eqs. (fTUjl. valid to zeroth order in e and e', 



A 



A +0 + V2A 



A 



+T) 1 



V2 

-^=A +0 + V2A 0v , 



07? 5 



(22) 



Measuring the magnitudes of the three amplitudes in each of the two triangles determines 
cosiprj and cos-i/y- Similar relations for charge-conjugate amplitudes determine if) v and i/y. 
To determine separately |A, o| and \A ( /)| would require measuring also CP asymmetries 



in these channels. In the absence of these asymmetry measurements, one may use charge- 
averaged rates alone to improve the upper bound ()19|). Maximizing (Aa — Aao)^- r?-r/ m 
(JTHJ) by varying if) o and Vyo, while keeping B + o and the upper bounds on B 0r){ i) fixed, we 



find that a maximum is obtained for \A 



\A 



0j7(')I : 



| (Act- Aa ] 



7T — TJ — 1J 



< 



where 



7-Q 



B 



+0 



B + , 



1 - r r , 



(23) 



16 



B +n B 



ry 



B+qBo v i 



Using world averaged values 



= (4.8 ± 0.6) x 10" 6 , B +rf 
and values in (j20J), we find at 90% CL 

|(A«- Aa ).- W | < 1.4° 



(4.2 ± 1.1) x 10~ b [3 



(24) 
(25) 

(26) 



This is only a slight improvement relative to (}2~T]) . 

The upper bounds (j2Tj) and (J2I3J) involve an uncertainty of about 30% from SU(3) breaking 
and small annihilation amplitudes which we have neglected. The bounds are seen to be 
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considerably lower than the estimate of the uncertainty, 5a ~ 5°, obtained in j3| using 
generalized factorization. These bounds may be tightened further by reducing errors in the 
relevant B + decay branching ratios, and in particular by improving the upper limits on 
BiB — > 7r°?7) and BiB — > 7r°r/). These experimental upper limits play also an important 
role in interpreting theoretically |4o| the measured deviation of the time- dependent CP 
asymmetry in B° — > rj'Ks from sin 2/3 sin Amt J4l| • This makes the case for their further 
improvement even stronger. 

Our analysis was based on the "ideal" mixing angle, 9o = sin _1 (— 1/3), which we used in 
(jHJ). Defining a general mixing angle, 9 = 9q + 5, one may show that for this case one must 
replace eo in (fT5|) and fl"6|) by e, 

[2 fl 
e = y -(ecos5 + e'sin<5) + y -(— esin<5 + e'cos<5) . (27) 

The small theoretical uncertainty in the value of 9 < 6°) implies a value for e within 
the uncertainty in eo given in ffT^j) . Furthermore, the terms e and e' in ()16|) are preserved by 
replacing 9q by 9, implying that the upper bounds (|2*Tj) and (j2HI) are unaffected by varying 
9. 



III. EFFECTS OF p - u MIXING IN B -> pp 

The processes B — > p % \n p 3 ^(^3^4) are quasi two-body decays involving four pions in 
the final state. To account for the p width, the two p mesons are defined by choosing 
suitable common ranges of invariant masses for the two-pion pairs, S12 = (pi + P2) 2 and 
S34 = (j>3 + P4) 2 . One uses the pion angular distributions in the p rest frames to project 
longitudinally polarized states which were shown to dominate B PP^%- Applying the 
isospin analysis to B — ► (pp)i ons proceeds identically to B — > 7T7T [ll, |2J in the limit of a 
vanishing p width. (In principle, the method applies separately to each transversity state.) 
The p width has the effect that two p mesons with different invariant masses, s±2 7^ S34, 
cannot be considered identical. Therefore Bose symmetry does not exclude a final 1 = 1 
state 0|, for which the amplitude is antisymmetric under su <-> S34. This amplitude does 
not interfere in the decay rate with the usual symmetric 1 = and 1 = 2 amplitudes. The 
effect of the 1 = 1 amplitude on the isospin analysis, of order (r p /m p ) 2 ~ 0.04, may be 
taken into account by including it in the fit. In principle, the effect may be eliminated by 
decreasing the width of the p band, however this would also decrease the statistics. 

In the following we will disregard this 1 = 1 term, which contributes also in the isospin 
symmetry limit, studying isosp in-breaking effects of the same order. To make our point, 
consider first the general invariant mass dependence of decay amplitudes for the three distinct 
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charged p states, 

^ + -(si 2 ,s 34 ) = A ( B ° - ^+n°) 12 (n\-) u ) = A(B° - p + p-)f c (s 12 )f c (s u ) , 

^4+0(^12, s u ) = A(B + -> (7T + 7r ) 12 (7r + 7r-)3 4 ) = -> p + p°)f c (s 12 )f n ( Su ) , 

Ax)(si2, s 34 ) = ^(#° -> (^ _ )i2(vr + 7r-) 3 4) = A(B° -> p p )f n (s 12 )f n (s u ) , (28) 

where / Cj „ are usually taken as Breit-Wigner factors. If isospin symmetry were exact, then 
fn( s ) — fc( s ), so that the two ratios, 

44^4 ^£^4, (29) 

would be independent of si 2 and S34 in the quasi two-body approximation. Any observed 
dependence of these ratios on the invariant masses would indicate either isospin-breaking, 
or dependence of A (B — > p % p>) on si 2 and S34. The latter possibility may be fitted exper- 
imentally by considering this dependence over the entire widths of the two p mesons [42]. 
We will study isospin-breaking in a narrow range of invariant masses defined by the narrow 
uj resonance, for which A(B — ► p l p*) may be assumed to be constant. Our purpose is to 
use the measured invariant mass dependence in (|2*Hj) as a tool for extracting the isospin 
symmetric B — ► pp amplitudes which obey a triangle relation similar to (UJ. 

Let us now study p — u mixing following a formalism developed in [43] . The physical p 
and uj fields are mixtures of an isovector field, pi, and an isoscalar fields, Uj, 

oj — uj + e 2 pi ■ 

The isospin-breaking parameters, ei j2 , are of order of a few percent. A precise knowledge of 
their magnitudes will not be needed (see appendix |X] for current experimental constraints), 
as they will be hidden in an isospin-breaking function to be introduced below. An expansion 
in 6i 2 will be carried out to first order in these parameters. 

Consider the transformation between the isospin basis and the physical basis for the 
scalar parts of the vector meson propagators. The mixed propagator in the isospin basis, 
= (pi^i)o, has poles at the p and u masses, and is conventionally written in the form 

D I puj (s)=Tl pw (s)D pp (s)DUs) ■ (31) 

The physical basis is defined by requiring that U puJ does not have poles. The scalar parts of 
the physical propagators can be approximated near the poles by Breit-Wigner forms, 

D ^ = s - K li mv r v ' v = ^- (32) 

The values of ei >2 are chosen such that the mixed propagator in the physical basis, D piJJ = 
{puj)o has no poles, 

D 1 ^ = + e 1 D uu - e 2 D pp . (33) 
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All three terms on the right-hand-side are of order ei i2 . The equalities, D vv (s) = Dyy(s), 
(V = p,w), hold to first order in ei )2 . For instance, the second and third terms in the 
relation, 

D pp = D pp + 2 £i D pu + , ( 34 ) 

are second order in ei j2 and will be neglected. 

To introduce isospin-breaking most generally, we take for neutral and charged p mesons 
independent p — > txix couplings, gi = g(pi — > 7i + n~), g c = g(p + — > 7r + 7r 3 ), and independent 
mass and width parameters entering D pp and D c pp . We neglect higher order effects in g(p + — > 
7r + 7T°) caused by 7r° — r] — rf mixing 



g{p + -> 7r + 7T°) = g(p + -> 7r + 7r 3 ) + eg(p + -> 7r + ??) + e'#(p + -> tt+t/) 



(35) 



because the two couplings multiplying e and e' violate G-parity and are thus further sup- 
pressed; for instance 0| 



g(p + -> 7T+77) 



#(p+ -> 7T+7r 3 ) 



m 
m 



Br(p + — ► 7r + 77) 



1/2 



< 0.055 (84%CL) . (36) 



5r(p+ 

In the presence of isospin-breaking u;j couples to two-pions with a coupling g(u;/ — ► nn) 
of order €ig(pi7nr). The decay B — > (7r7r) X then proceeds either through pj or through uj. 
Working to first order in isospin-breaking, these two contributions enter through a linear 
combination of the two propagators, both of order e*, 



D 1 



7T7T 



g(pi 



-d 1 



1T1T 



(37) 



Thus, one finds expressions for B decay amplitudes into four pions, including terms which 
are first order in e,-: 



— ( s 12> S34) 
A +0 (Si2, s 34 ) 



g c gi 



P + P')D c pp (s 12 )D c pp (s u ) 



(38) 
(39) 



+ 7=2 A{B ° 



A(B + -> p + Pl )D pp {su) 
+ A(5+ - p + ^)^( S34 )]z^ p ( Sl2 ) , 
^oo(si2,s 34 ) = #f ^(5° -> pipi)D pp (s 12 )D pp (s u ) 

Piuji)(d pu (si2)D pp (s 3 4) + (S12 ^ s 34 )) . (40) 

An implicit angular dependence in ()38j) - (j40jl . corresponding to given polarization states 0], 
is independent of si 2 and s 34 . 

In the isospin symmetry limit, D c pp = D pp , D puj = 0. Isospin-breaking is given by devia- 
tions from these equalities and, for the case of p — u mixing, is parametrized most generally 
by Eqs. ()38j) - (j40j) . Each term in a given row has a distinct dependence on si 2 and s 34 , char- 
acterized near the p and u poles by (J3TJ), and (jHTJ). Taking g c jg 1 = 1.005 ± 0.010 |4J, 
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the invariant mass distributions of the three processes permit in principle a determination 
of the three magnitudes, \A(B° — ► p + p~)|, \A(B + — > p + p/| and \A(B° — > pipi)\, forming the 
isospin triangle, 

-> p+p~) + V2A(B° -> - ^(5+ -> p+p 7 ) = . (41) 

Once this triangle and its charge-conjugate are formed, one uses a phase relation for A^B^ — > 
p ± pi) analogous to (JBJ) and the CP asymmetry in -B°(t) — > p + p~ to determine a. This then 
provides a way of extracting a free of effects from p — u mixing. Electroweak penguin 
contributions are treated as in B — > 7r7r, Eq. (JZj). 

The extraction of the pure isospin amplitudes \A(B — > pp)| may be facilitated by using 
information from direct measurements of A(B + — > p + cu) and — > p°c<j) entering the 
isosp in-breaking terms in ([3~S j) -PD j) . Also, the isospin-breaking function D^s) is the same 
as the one fitted to the pion form factor |46j . Denoting 

1 

[s-m 2 p + irripTp] [s - ml + im w T u 

t2 

,2 TT' 1^2- 



D P M) = rV00 r „ Z3 I TZ i [„ > ( 42 ) 

the fit yields il ptlJ (m^) = —3500 ± 300 MeV 2 , involving a possible small imaginary part 



compatible with zero. The slope at s — , II' (to 2 ) = 0.03 ± 0.04, is consistent with zero. 
The exact s dependence of U puJ is unimportant because its contribution is dominated by the 
narrow uo width. 

At the uo mass this gives 

\Dpu{rnl)\ 
\D pP {ml)\ 

Using the experimental values 



0.53 ±0.05. (43) 



B{B + -> p+p°) = (26.4±li) x 10" 6 , B{B + -> p+c^) = (12.6±|;|) x 10~ 6 , (44) 

and neglecting possible CP asymmetries in these processes, leads to 

\A(B + ^ p + uo)D p Jm 2 )\ . . 

j-^— 4- n ; = 0.36 ± 0.08 . 45 

While this ratio becomes 0.02, typical for isospin-breaking, when weighed by the uo and p 
widths, it has a large effect at the u> mass. 

In order to demonstrate the effect of p — uo mixing on the 7r + n~ invariant mass distri- 
bution applying Eq. (|39j) . one must use some information about the relative magnitudes 
and relative phases of the amplitudes for B^ 1 — > p ± p° and B + — > p^. While the former 
amplitudes are pure "tree" (we neglect very small EWP contributions), involving a single 
CKM phase tKxg(V* b V ud ) = 7, the latter involve also penguin contributions with weak phase 
Arg(V* b V cd ) = ir in the c-convention 0], 

V2A(B + -> p+p°) = t + c , ^,4(5+ -> p+cj) = t + c + 2p + 2s . (46) 
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FIG. 1: Invariant mass distributions for 7r + 7r~ in B — > p it + tt~ demonstrating p — uj mixing, 
using hadronic parameters as given in the text. Dashed (blue) line represents B + decays; solid 
(red) line represents B~ decays; thick (black) line describes a case neglecting p — u> mixing. 



For our purpose, the terms t, c, p and s represent SU(3) amplitudes for longitudinally polar- 
ized vector mesons, similar to those defined in Sec. II for two pseudoscalars. The amplitude 
s is OZI-suppressed and is expected to be negligible. (A similar amplitude in decays to a 
vector meson and a pseudoscalar meson dominates B + —>■ (f)7T + [50].) The ratio |p|/|t + c| 
is small, about 0.1, as can be inferred from the small branching ratio of B° — > p°p° 6||, or 
(by flavor SU(3)) from the measured longitudinal branching fraction of B — >• K*(p 33]. 

For illustration, we use 2|s + p|/|t + c| = 0.2, 7 = 57°, choosing the strong phase difference 
between s + p and t + c to be zero, so that the ratio B(B + — ► p + u)/B(B + — ► p + p°) = 0.82 
is smaller than one as implied by experiment, Eq. (|44|). (Small transverse contributions are 
neglected.) These parameters determine relative magnitudes and relative phases between 
A(B ± — > p ± oj) and A{B ± —>■ p ± p°). The resulting effect on the ti + tt~ invariant mass 
distribution is shown in Fig. 1, separately for B + and B~ decays. The vertical axis gives 
number of events in arbitrary units. The prominent peak, followed by a dip, is characteristic 
of p — <jj mixing |5lJ , and does not depend strongly on the choice of hadronic parameters as 
long as |p + s|/|t + c| <C 1. Note that a small CP asymmetry is expected also in case that 
the strong phase difference between s + p and t + c vanishes, because of the two different 
shapes of the p and u resonances [3, [3, E3] • 

In reality, limited statistics and particularly the current absence of a positive signal for 
B° — ► p°p° would forbid carrying out the complete program leading to a construction of 
the pure isospin triangle (|41|). However, using the given invariant mass dependence of the 
isosp in-breaking terms in ()38j) - (j40p . one may constrain these terms and eliminate them in 
certain cases. As noted above, a first place to look for these terms would be the 7r + 7r~ 
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invariant mass distribution in B + — > p + n + n~ near the u mass, which should be fitted to a 
sum of the p + p° and p + u terms in (J39j) as plotted in Fig. 1. 

In the relativistic Breit-Wigner form (|32|) we assumed no s dependence in the width 
IV. This assumption has only a very slight effect on the isospin-breaking terms which are 
dominant at the narrow u peak. Non-resonant contributions could also affect the invariant 
mass dependence. One hopes to minimize these contributions by fitting the pion angular 
distributions to those describing longitudinally polarized vector meson states. Potential 
interference between the p resonance and the wide p'(1450) resonance must also be taken 
care of. 



IV. MESON MIXING AND OTHER ISOSPIN-BREAKING EFFECTS IN B -> pir 

The isospin method for extracting a in B — > pn is based on a time- dependent Dalitz plot 
analysis of B° — > 7r + 7r~7r° [5;], using information provided by isospin symmetry We will 
explain now the essence of the method. Consider the contributions of p + ,p~ and p° to the 
amplitude of B° — > 7r + 7r~7r°, 

A + _ = A(B° -> tt+tt^tt ) = A + D pp {s + ) cos 6 + + A^D PP (s_ ) cos 9^ +A D pp (s ) cos9 , (47) 

where subscripts of the amplitudes A iy the invariant masses Sj, and the helicity angles 6i, 
denote the charge of the p, 

A + = A(B° -> p+vT) , A. = A(B° -> p-7r+) , A = -> p°vr ) . (48) 

The function D pp is given near the p pole by a Breit-Wigner form (|32p. Corresponding 
amplitudes for B° are denoted by A. The time-dependent decay rate for an initially B° is 
given by [^l 

r(jB o _> ,r+7r-7r (t)) oc (|A + _ | 2 + |A_ +0 | 2 ) + (|A + _ | 2 - |A_ +0 | 2 ) cos(Amt) 

- 21m (e- 2i/3 A_ +0 A;_ ) sin(Amt) . (49) 

The interference of the three p resonances in the time-dependent and invariant mass- 
dependent decay rate permits a determination of the magnitudes of the three amplitudes 
A + , A_, A and their charge conjugates, as well as the relative phases between these six 
amplitudes. This amounts to eleven independent observables, obtained from a set of twenty 
seven mutually dependent measurables in PH)) [5J. It is convenient to rescale the A\ and Ai 
amplitudes by phases e t/3 and e~ zlS , respectively, defining Ai = exp(i{3)Ai, At = exp(— i0)A\, 
such that the coefficient of the sin(Amt) term in measures directly the relative phases 
between Ai and A4. 

Working in the t-convention 0|, each amplitude Ai consists of a "tree" contribution 
proportional to V* b V u d and a "penguin" term involving V t * b Vtd- Unitarity of the CKM matrix 
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is used to absorb matrix elements of penguin operators proportional to V* b V u d in the tree 
amplitude, 

At, = e- 4a T ±i0 + P ±t0 , A ±fl = e +lQ T ±i0 + P±, , (50) 

where T ± and P± t0 include strong phases. The measured left-hand sides provide eleven 
equations for twelve unknown parameters, consisting of the weak phase a, six magnitudes of 
tree and penguin amplitudes, and five relative strong phases between these amplitudes. An 
additional complex relation, leading to a total of thirteen equations for the twelve param- 
eters, is provided by isospin symmetry. Neglecting electroweak penguin contributions and 
isopin breaking effects, the pure AI = 1/2 penguin terms vanish in the 1 = 2 amplitude, 
implying 

P_ + P+ + 2P = . (51) 

The over-constrained set of equations ([50|h ()51|) allows for the extraction of a. The sensitivity 
to a can be seen explicitly by considering the AI = 3/2, 1 = 2 amplitude A2, 

A 2 = A+ + A- + 2A = Te~ ia , A 2 = A+ + A- + 2A = Te ia , (52) 

with T = T + + T_ + 2T . The angle a is fixed by the relative phase between the two sums 
of amplitudes, which are determined up to an overall common phase by the time- dependent 
Dalitz plot. 

The relations (|51|) and ()52)1 are violated both by EWP contributions and by isospin- 
breaking corrections caused, for instance, by 7r° — r] — rj mixing and p — lo mixing |l7j. 
The former contributions may be included model-independently as in the case of B — >• tttt 
and B — ► pp, using the proportionality of the AI = 3/2 current- current operator and the 
dominant AI = 3/2 EWP operator in the effective weak Hamiltonian The relation 

(foTj) is modified to 

P- + P+ + 2P = P EW , (53) 

where Pew can be obtained from 

A 2 = Te- + Pew , ^ = -\( ^) ^ = + 0.013^^i . (54) 

T 2 \ ci + c 2 / smfj 

This implies the same shift in a as in B —>■ tttt, Aoewp = (1.5 ±0.3)°, if only the information 
from the two measured sums of amplitudes A2 and A2 is used to extract a. The shift in a 
obtained from fitting the entire set (J3UJ), and (fSlj) may be slightly different, because of 
the additional a dependence of this over-constrained system of equations. 

Isospin-breaking in tree amplitudes does not affect the extracted value of a in B — » 
pn, which is based on Eqs. (IHUj) . and because no isospin relation between tree 

amplitudes is needed. This is obvious when a is extracted from (1521) and is true in general. 
Penguin aI nplitnda S in £? — * p7T a r e known to be small H ^Therefore, other SO ntce S 
of isospin-breaking are expected to lead to corrections in a smaller than Aoewp; which is 
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related to the 1 = 2 tree amplitude through The evaluation of corrections caused by 

7T° — i] — r)' mixing, that we give next, supports this expectation. 

The mixing of 7r°, rj and rj' introduces additional correction terms in ([52)1. 

A+ + A- + 2A = Te~ ia + P EW + eP pv + e'P m , , (55) 

where P prj C) are penguin amplitudes in B° — > p°rf\ We expect that the correction to the 
extracted value of a is smaller in B — > p7r than in i? — > 7r7r because the p eng uin-to-tree ratio 
is smaller in the first case. In terms of SU(3) amplitudes defined in |50l l55j]. one has 



P m = —(-Pp-Pv ~ s v ) 



\T\ = \t P + t V + C P + Cy| 
1 

P PV ' = ~^f Pp +Pv + As v) ■ (56) 
A global SU(3) fit to available data of charmless B decays to a pseudoscalar and a vector 



meson 



55j has shown that tp and ty add up constructively, while cp and cy are smaller. 



P a nu by cH 

Also, \r>p/tp\ ~ \pv/tv\ ~ 0.2 0, and pv — —pp, while sy is smaller. (A best fit 
gives \pv/pp\ = 1.15 ± 0.07, arg(p y /p P ) = (182 ± 18)° and sy/py = 0.16±g;g|.) All this 
implies that the effect of the terms in (|55J) involving e and e' is very small. Taking 

\T\>\t v \, \P m \~^\s v \<^\p v \, [P^l ~ A| Sv | < 2|| pv | , ^ = 0.2,(57) 

and using e = 0.017 ± 0.003, e' = 0.004 ± 0.001, we find the following upper bound on the 
uncertainty in a caused by neglecting the terms in (|55|): 

| Aa,_„, | = |ePpr? ^ /Pp? ' ,| < 0.024e + 0.069e' < 0.1° . (58) 



In case that the sum of pp and py in P prj C) does not cancel completely |55[, the bound could 
be a factor two larger. In any event, this uncertainty is much smaller than A«ewp; the shift 
caused by EWP amplitudes. 

Finally, we discuss the effect of p — u mixing treating it as in Sec. III. Neglecting isospin- 
breaking in g(p — > mi), we conclude that the third term in f)47|) must be replaced by 

A (B° -> p\°)D pp (s ) - A(B° -> ^^(so) + A(B° ^ Wjtt )^^) , (59) 

while the angular dependence remains unchanged. That is, the effect of p — u mixing 
may be included in the time-dependent Dalitz plot analysis by adding the second term 
in ()59|) . The isospin-breaking function D pul (s ), defined in (J37j) and given in (|4*2^) . has a 
double pole structure with a narrow peak at the u mass. As discussed in the previous 
section, D puJ (s ) is measured by studying the pion form factor, while A(B° — > Ldfii ) can be 
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obtained from B — > Air. The narrow peak at the 00 mass distinguishes clearly this isospin- 
breaking correction from other potential non-resonant or wide resonance contributions to 
B ^ pn [llEil- 



The size of the effect of p — u mixing may be estimated by considering the two processes, 
B° — > p°7T and B° — > um°, occurring in ()59|1. The charge-averaged branching ratio of 
B — > p°7T measured by Belle is somewhat larger than an upper limit reported by BaBar, 
reporting also the currently strongest upper bound on B° — > um°, 



B(B° - p ^ 



(5.1 ± 1.6 ±0.9) x 10~ 6 , Belle [gg] , 
< 2.9 x 10~ 6 , BaBar [5^] , 
£(5° -> cj7t°) < 1.2 x 10~ 6 , BaBar . (60) 



These values and permit a relatively sizable contribution from the isospin-breaking term 
A(B° — > ujiir )D f)UJ {.so) at the pole, so = m^, if — > ujtt°) is not much below its current 
upper limit |55| . 



V. OTHER SOURCES OF ISOSPIN VIOLATION 



In this work we have focused primarily on isospin-breaking effects in B — ► ttti, B — > pp 
and B — > p7r, originating in the mixing of neutral isospin triplet states (7r°,p°) with isospin 
singlet states {rf\uS). We also iterated the effects of higher order electroweak penguin 
operators. Two implicit assumptions were made in our analysis: 

• Reduced matrix elements of operators in the effective Hamiltonian, between initial B° 
and B + states and final states involving 7r 3 and 7r + , were assumed to obey exact SU(2) 
relations. 

• AI = 5/2 corrections were assumed to vanish. 

Relaxing these assumptions in B — > titt and B — > pp introduces isospin- violating corrections 
in a which may be comparable to those discussed in sections II and III. AI = 5/2 operators 
in B — > 7T7T jl3 or in B — >• pp may be induced by an insertion of the d — u mass difference 
AI = 1 operator or by electromagnetic corrections. This would violate the closure of the 
isospin triangles for B and B amplitudes. Any of these other isospin-breaking corrections in 
B — » p7r is expected to be negligible, however, because in these processes isospin breaking 
can only affect the relation among suppressed penguin amplitudes. 

To see how these other effects enter a specific calculation, let us use the result of a 
Soft Collinear Effective Theory approach to factorization in B — » MiM 2 , where M 12 are 
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(61) 



pseudoscalars or vector mesons. The result, at leading order in A/mg, is jl^ |: 
A = ^f^M f dudzT^u, z)(f M > {z)<j> Ml {u) 

+ f Ml C BM2 f T K (u)^(u)} + {l - 2} + \i»A^ M > . 

Here Tjj(w) and T^u) are hard kernels which may be expanded in ^(wb), while ( BM ', 
(f M {u) and the light cone meson wave function <fi M (u) are nonperturbative parameters. 
The amplitude A^' M2 denotes a possible long distance charming penguin contribution. 

We have quantified isospin-breaking caused by final states which do not coincide with 
isospin eigenstates. The remaining isospin violation is encoded in ( BM , ( BM (z), fM<t> M {u), 
and A^ 1,Mi , where are now isospin eigenstates. Generically, the corrections are expected 
to be of order (m u — m^/^-QCD ~ « ~ 0(1%), namely of the same magnitude as the 
corrections caused by n — 77 — rf mixing, Eq. (j2HI), and by EWP contributions, Eq. (JZj). 
Isospin-breaking in the hard kernels Tij(u),T i( ^(u) occurs through additional 1/mg power- 
suppressed operators and may be safely neglected. At this order, no isospin violation is 
caused by final state rescattering in the first two terms in the amplitude ()61|) . because these 
terms factorize to all orders in as and to first order in A/mg. 



VI. CONCLUSIONS 



The extraction of the weak phase a = 02 by application of isospin symmetry to B — > ttty, 
B — ► pir and B — ► pp is modified through 7r° — 77 — 77' mixing and p — uo mixing. We 
have studied these effects in a model-independent manner, discussing also other effects of 
isospin-breaking in these processes. Our main results are the following: 

• Isospin-breaking corrections in a related to tt° — 77 — rj' mixing were bounded using 
flavor SU(3), and were found to be smaller than 1.4° in B — > ttti and much smaller in 
B — > pir. 

• The effects of p — uj mixing in B — > pp and B — > pir were studied as a function of the 
two-pion invariant mass in terms of a quantity measured in the pion form factor. Given 
the invariant mass dependence characterizing p — u mixing, which involves a peak at 
s = ni^, we propose a way for measuring and constraining these effects experimentally. 
Eventually, with sufficient statistics, this procedure may eliminate the mixing effect 
altogether. 

• In B — > p7r, any kind of isospin-breaking in tree amplitudes does not affect the measure- 
ment of a through a time-dependent Dalitz plot analysis. Since penguin amplitudes 
are suppressed, the resulting uncertainty in a from isospin violation is expected to be 
smaller than one degree (excluding contributions from EWP operators). 
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• The proportionality of a AI = 3/2 current- current operator and a corresponding 
dominant electroweak penguin operator in the effective Hamiltonian implies a shift, 
Aa EWP = (1.5 ± 0.3)°, common to B — > nn and B — ► pp. The same shift would apply 
also to B — > pit if only the sums of amplitudes (|52j) were used. In a completely general 
extraction of a from the time-dependent Dalitz plot fit, the shift may be slightly 
different but can be obtained model-independently. 

A brief summary of our conclusions is therefore: (1) Isospin-breaking introduces a much 
smaller uncertainty in the value of a extracted from B — > tttt than thought before, of order 1°. 

(2) Effects of p — u mixing in B — > pp can be studied by fits to invariant mass distributions. 

(3) The largest shift in a in B — > pn, caused by electroweak penguin amplitudes, can be 
included model-independently, and is about 1° as in B — > tttt and B — ► pp. 
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APPENDIX A: EXPERIMENTAL CONSTRAINTS ON e 1>2 

Let us comment briefly on the values of the isospin-breakin g pa rameters, e^, which are 
constrained by fitting D puJ (s) in (|4^j) to the pion form factor [46|. Requiring that D puj (s) 
does not have poles at m 2 — im p T p and — im^Y^ implies 

= iWmg, - im^T^) ^ = U^jml - im p T p ) 

Using the relation 

fi^ (s) = Ii pu) (s) + 9 Y* ~" ™\ (s -m 2 + im p T p ) , ( A2) 

one may express e« in terms of the measurable function n pu; (s): 



ei 



Rpuirnl - imj: u ) gfa tttt) 



ml-m? p + i(m p T p - m u r w ) g(pi -> tttt) 

ti pu) {m 2 - im p T p ) 
m^-mi + i{m p T p - m u T u ) 
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The first term in e\ is constrained experimentally 46] (see also discussion below 

' [/ '-' ( "' ; ~ "" J ~' ) - (-0.003 ± 0.002) + i(0.032 ± 0.003) . (A4) 



Barring a possible weak dependence of Il pw on s, this term is equal to e-i- The term 
g{uj — > irir)/g(p — » 7r7r) in ei is poorly constrained experimentally, but is expected to be 
of the same order. The smallness of these parameters justifies neglecting terms of order ef 2 - 
Note that the method presented in Sec. Ill for studying isospin-breaking in B — > pp depends 
on the function Yl puJ (s) and not separately on its two components given in (|A2|) . 
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